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The Yukawa-Higgs/Ginsparg- Wilson- fermion construction of chiral lattice gauge theories described 
in hep-lat/0605003 uses exact lattice chirality to decouple the massless chiral fermions from a 
mirror sector, whose strong dynamics is conjectured to give cutoff-scale mass to the mirror fermions 
without breaking the chiral gauge symmetry. In this paper, we study the mirror sector dynamics of 
a two-dimensional chiral gauge theory in the limit of strong Yukawa and vanishing gauge couplings, 
in which case it reduces to an XY model coupled to Ginsparg- Wilson fermions. For the mirror 
fermions to acquire cutoff-scale mass it is believed to be important that the XY model remain 
in its "high temperature" phase, where there is no algebraic ordering — a conjecture supported by 
the results of our work. We use analytic and Monte-Carlo methods with dynamical fermions to 
study the scalar and fermion susceptibilities, and the mirror fermion spectrum. Our results provide 
convincing evidence that the strong dynamics does not "break" the chiral symmetry (more precisely, 
that the mirror fermions do not induce algebraic ordering in two-dimensions), and that the mirror 
fermions decouple from the infrared physics. 
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1 Introduction and summary 



1.1 Motivation 

The study of strong-coupling chiral gauge dynamics is an outstanding problem of great interest, 
both on its own and for its possible relevance to phenomenology. Whereas the standard model 
of elementary particle physics is a weakly coupled chiral gauge theory, additional strong chiral 
gauge dynamics at (multi-) TeV scales may be responsible for breaking the electroweak symmetry 
and fermion mass generation. Various tools are currently available for the study of the strong- 
coupling behavior of chiral gauge theories. For instance, one has 't Hooft's anomaly matching and 
most attractive channel arguments. These are complemented by the "power of holomorphy" in 
supersymmetric theories. Scaling arguments and effective NJL-like models, both using results from 
QCD as a stepping stone, have also been employed extensively. (Large- iV expansions, including the 
recently considered gravity duals in the AdS/CFT (AdS/QCD) framework, do not usefully apply 
to chiral gauge theories.) None of these approaches represents a "first principles" method, with an 
accuracy that can be systematically improved. Thus, the space-time lattice regularization remains, 
to this day, the only way to advance our limited knowledge of strong chiral gauge dynamics. 

During the past two decades, since the work [1] of Ginsparg and Wilson (GW), there has been 
significant progress in understanding chiral symmetries on the lattice [2-6] (further references are 
given in the reviews [7,8], while [9,10] contain more recent work). However, the lattice formulation 
of chiral gauge theories is still a largely unsolved problem, although there were some advances 
and promising two-dimensional numerical results using the original overlap Weyl determinant [11]. 
The fact that a general formulation is lacking, along with the natural expectation that the lattice 
definition of a chiral gauge theory is not unique, indicates that the exploration of the various 
existing formulations, of their interconnections, as well as of novel constructions of chiral lattice 
gauge theories are both warranted and worthwhile. 

This paper is devoted to a study of a recently proposed construction of lattice chiral gauge 
theories [12]. It combines the GW exact lattice chiral symmetry with earlier ideas of strong coupling 
Higgs- Yukawa dynamics on the lattice, see [13-17] and references therein. The essence of the 
proposal is that the exact lattice chiral symmetry may allow an old dream to be realized: that 
the mirror partners of the chiral massless fermions can be decoupled, without breaking the gauge 
symmetry. The decoupling of the mirrors is possible in the strong- Yukawa symmetric phases of 
lattice fermion-Higgs models, where the fermions have cutoff-scale mass without breaking the chiral 
symmetry. 

The mechanism of mass generation without chiral symmetry breaking operative in strongly- 
coupled lattice theories has been known for some time^l at strong Yukawa coupling, the charged 
mirror fermions form neutral bound states with the charged scalars; these bound states can pair 
up with other neutral mirror fermions, composite or elementary, to obtain mass without breaking 
the chiral symmetry. 

The construction of [12] gives an explicit definition of a local lattice action and measure. The 
global symmetry Ward identities, exact or anomalous, are precisely as in the target continuum 
theory, which sets the proposal apart from earlier constructions aiming to decouple the mirror 

In the two-dimensional context this mechanism was elucidated by Witten [18] via bosonization. In the case of 
four-dimensional strong four-Fermi interactions on the lattice, it was part of the proposal of Eichten and Preskill [13] 
to formulate chiral gauge theories on the lattice. In the context of strong- Yukawa lattice models it was discussed 
in [14-17]. 
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fermions. We believe that these desirable features are alone sufficient to motivate further study of 
the proposal. 

The lattice model of [12] will give rise to an unbroken chiral gauge theory only if there exists a 
strong- Yukawa phase with an unbroken gauge symmetry, massless chiral charged fermion spectrum, 
and massive mirror fermions. Because the construction makes use of the somewhat complicated 
exact lattice chiral symmetry, implemented via the Neuberger-Dirac operator [4], a strictly analytic 
approach — e.g., a strong-coupling expansion as used in [13-16] — to establishing the existence of 
the desired strong-coupling phase appears out of reach. The purpose of this paper is to provide 
evidence, using a semi-analytic approach that computes the leading terms of the strong-coupling 
Yukawa expansion by Monte Carlo estimation, for the existence of the strong-coupling symmetric 
phase with a massive mirror fermion spectrum. 

1.2 Organization and summary of main results 

In Section [21 we begin with a review of the proposal and of the arguments leading us to expect that 
an unbroken gauge theory with a chiral spectrum of charged fermions emerges in the infraredH 
We present the construction on the example of a two-dimensional chiral U(l) theory, the well- 
known "345" theory. While the proposed lattice action can also be written in four dimensions, 
the simplicity of the analytical analysis, and especially, the relative ease of the numerical analysis, 
restricts the present study of the dynamics to the two-dimensional case. 

Of course, due to the special properties of two-dimensional theories, the exact U(l) symmetry 
cannot be spontaneously broken [19]. However, the gauge boson would acquire a large mass (relative 
to the gap induced by the Schwinger mechanism) in the quasi-ordered phase, which is analogous to 
spontaneous symmetry breaking in the four-dimensional case. It is for this reason that here, and in 
several places below, we shall, with some abuse of language, use the term "broken" phase to refer 
to the low-temperature phase of the two-dimensional XY model, instead of the more appropriate 
algebraically- or quasi-ordered phase. Similarly, we will sometimes refer to the high-temperature 
disordered phase as the "symmetric" or "unbroken" phase and will take the liberty to omit the 
quotation marks. 

In Section [3l we formulate a "toy" version of the model, simple enough to be subjected to ex- 
tensive analytical and numerical tests. The main focus of this paper is the Higgs- Yukawa dynamics 
of this "toy" theory at strong Yukawa coupling and vanishing gauge coupling — the relevant limit to 
approach the continuum theory. The importance of this analysis is that it addresses precisely the 
question: is there a symmetric phase with a massless chiral spectrum of charged fermions in this 
model? A positive answer to this question is crucial for the viability of the proposal; a negative 
one would imply that the proposal should be abandoned. 

We begin the analysis in Section 13.2^ where we show that the partition function of our model 
factorizes into a product of "light" and "mirror" partition functions. We work out the splitting of the 
fermion measure into a product of positive and negative GW-chirality measures, using Neuberger's 
Dirac operator, and the resulting factorization of the partition function; we give the technical details 
pertinent to this analysis in Appendix [A] The factorization of the partition function is a unique 
feature of the GW-fermion Higgs- Yukawa theory that makes it distinct from other Higgs- Yukawa 
approaches, notably the "waveguide" models [20-22]. 

2 See also Section IV in [12]. 
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The exact light/mirror split of the partition function at g = allows us to further concentrate 
entirely on the study of the dynamics of the mirror fermion-Higgs theory. The g = mirror fermion- 
Higgs theory of the "toy" model is equivalent to the two-dimensional XY model (see, e.g., [23]) 
in the high-temperature disordered phase, k < K c , deformed by the addition of fermions strongly 
coupled to the spins. 

We next ask, in Section 13.31 whether a strong coupling symmetric phase of the mirror-fermion 
Higgs- Yukawa theory exists. We begin by noting that at strong Yukawa coupling the mirror fermion 
determinant is positive for arbitrary Higgs field background, provided the ratio of the Majorana to 
Dirac Yukawa couplings h > 1. The positivity of the fermion determinant permits us to numerically 
study the effect of mirror fermion "loops" on the Higgs field dynamics in the symmetric phase (i.e., 
at small Higgs kinetic term k)H 

We demonstrate the existence of a symmetric phase at strong Yukawa coupling by studying 
the scalar susceptibility in Section 13.3. Q the Binder cumulant — a quantity [24] probing higher- 
order correlations — in Section 13.3.21 the vortex density in Section 13.3.31 and the fermion bilinear 
susceptibilities in Section 13.3.41 We use numerical simulations to investigate the behavior of the 
various susceptibilities for different values of k and h on lattices of size N 2 = 4 2 , 8 2 , 16 2 . The 
Monte Carlo techniques we use are described in Appendix O The combined study of the scalar 
and fermion bilinear order parameters and their scaling with N presents ample evidence that the 
theory is in the symmetric phase at strong Yukawa coupling to the mirror fermions, provided h > 1. 

For h < 1, still at small k, we find evidence for criticality as appropriate in the low-temperature 
phase of the XY-model, where the susceptibility increases with the volume of the system; con- 
currently we find that the mirror fermion spectrum contains massless modes and that the vortex 
density decreases rapidly as h approaches unity from above, as in the Berezinski-Kosterlitz-Thouless 
transition to the low-temperature phase. 

In summary, our results of Section 13.31 show that the strong coupling symmetric phase exists 
and thus the U(l) symmetry — the global part of the gauge group — is unbroken. This result also 
implies that no fine-tuning is required in order to keep the theory in the symmetric phase (after 
gauging — keeping the gauge boson perturbatively massles^l) at strong Yukawa coupling. This is in 
contrast with a naive perturbative argument which would, at large Yukawa coupling, imply a large 
backreaction of the mirror fermions on the Higgs dynamics, thus requiring fine-tuning of potentially 
many terms to keep the model in the symmetric phase. The small backreaction of the fermions on 
the Higgs dynamics at large Yukawa is, however, consistent with the qualitative strong-coupling 
arguments of Section [2] and [12]. 

The next question we address is the mirror fermion spectrum in the symmetric phase. In Section 
13.41 we study the spectrum of both charged and neutral mirror fermions and show that they are 
all massive in the strong Yukawa symmetric phase, for h greater than unity. 

We note that, to the best of our knowledge, this is the first strong- Yukawa/Higgs construction 
within the framework of the "waveguide-like" (or Wilson- Yukawa) approaches to lattice chiral gauge 
theories that successfully yields a chiral massless spectrum of fermions at g = 0. Previous Yukawa 
approaches have failed to yield a chiral light spectrum of charged fermions already at the quenched, 
g = 0, level, see [21,22] and references therein. (The only exception that we know of is the quenched 
study [25] of a theory that is in some respects similar to our toy model, in the framework of the "two- 

3 Actually, "loops" is a bit of a misnomer since our large Yukawa expansion is just the opposite of the usual loop 
expansion. Our results are in that sense nonperturbative. 
4 Setting aside the effects of the Schwinger mechanism. 
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cutoff' construction of [26]. The advantages of the model that we study here is that, unlike [25], it 
exactly preserves the symmetries of the target theory and the mirror fermions exactly decouple.) In 
our construction, the crucial ingredient that leads to a chiral spectrum is the decoupling, due to the 
exact lattice chirality, of the light chiral fermions from the strong Higgs- Yukawa dynamics of the 
mirror sector. It is the main feature that distinguishes our models from earlier "Wilson- Yukawa" 
or "waveguide" constructions. 

1.3 Conclusions and outlook 

Our main results are: 

1. We have shown that exact lattice chirality can be used to decouple the massless chiral fermions 
from a mirror sector. The light-mirror split of the fermion partition function is given in 
Appendix [A] 

2. We have given numerical evidence that the strong Yukawa dynamics of the mirror sector gives 
cutoff-scale mass to the mirror fermions, without breaking the gauge symmetry (i.e. without 
inducing algebraic ordering in two dimensions). 

3. We have found that, at strong Yukawa coupling, the main effect of the mirror fermions on the 
XY model to which they are coupled is to renormalize the hopping parameter n to smaller 
values — deeper into the "high temperature" phase — provided h > 1. For instance, in the XY 
model coupled to fermions, the vortex density is higher, and the susceptibility is lower, than 
in corresponding pure XY model. This conclusion is further supported by the step function 
centered on h = 1 that develops in the Binder cumulant, extrapolated to the N 3> 1 limit. 

Our results show that our proposed formulation of chiral lattice gauge theories satisfies a first check 
on its viability: the strong Yukawa dynamics produces heavy mirrors without breaking the gauge 
symmetry. 

Clearly, the results of this paper encourage further inquiry: 

• The toy model studied here gives rise to a would-be anomalous massless fermion spectrum in 
a trivial (U = 1) gauge background. The issue of gauge anomalies in models where non-gauge 
dynamics in the mirror sector is introduced to decouple the mirrors was recently studied 
in [27]. As explained there, in the anomalous case a result for the mirror spectrum at U = 1 
can not be used to infer the spectrum at U ^ 1, since for an anomalous fermion content the 
mirror fermion partition function is not a smooth function exactly at U = 1 and the split 
of the partition function into "light" and "mirror" is not consistent in this case (this is how 
the non-gauge mirror sector dynamics introduced to decouple mirrors is required to obey the 
anomaly free condition). Clearly, the next important step is to study the mirror dynamics 
in an anomaly-free model, like the 345 model. See also the discussion in the Addendum of 
ref. [27]. 

• In the anomaly free case, we have reason to believe that our results will hold for nonvanishing 
gauge coupling, since the mirror- fermion/ X Y-model physics is at the ultraviolet scale, where 
the gauge coupling is weak. As shown in [27], the mirror partition function is a smooth 
function of the gauge background in the anomaly free case and a result showing that the 
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mirror fermions decouple at U = 1 should persist for small nontrivial gauge background, e.g. 
in perturbation theory with respect to the gauge coupling. 

• A study analogous to the one of this paper in the four-dimensional case is both feasible and 
desirable. We expect that many of the details will be different — this is clear already from the 
fact that in four dimensions the need for neutral massless "spectator" fermions, dictated by 
two-dimensional Lorentz invariance, does not arise. We note in this regard the study of [28], 
showing, within an analytic strong-coupling expansion the existence of a strong-coupling 
symmetric phase in a four-dimensional Yukawa-Higgs model with GW fermions (this result 
is backed-up by Monte Carlo simulations [29]). 

2 The proposal 

In this section, we review the proposal of [12]. It has many desirable features: 

• It is a full lattice proposal for a local action and measure for a chiral gauge theory. 

• The realization of both anomalous and anomaly-free global symmetries is exactly as in the 
target continuum theory. 

• Arguments in support of the conjecture that at strong Yukawa coupling the theory is in a 
symmetric phase with a chiral massless spectrum were given in ref. [12]. 

We believe that the three points above warrant the further study of the proposal. In this paper, we 
investigate the properties of this "Yukawa-Higgs-GW-fermion" lattice model in more detail. First, 
we repeat the formulation of the model [12]. Then, we subject the third point above to a detailed 
analysis. 

2.1 The continuum "345" model and its symmetries 

We will present the proposal on the example of a two-dimensional U(l) chiral gauge theory. The 
condition for a two-dimensional U(l) gauge theory to be free of gauge anomalies is that the sum of 
the squares of the charges of the left and right handed modes must cancel ^ Qi i e ft~^2j 1j right = ®- 
An example of such a theory is the "345" theory where there are left-handed fermions of charge 3 
and 4 and right-handed fermions of charge 5, so the chiral fermion spectrum is simply denoted as 

3 , 4_, 5+. The model is asymptotically free and can be solved via bosonization (its spectrum was 
found in [30]). 

The "345" U(l) theory has two anomaly free global symmetries, the "133" symmetry (in a 
notation where 3_ has charge 1, and 4_ and 5+ — charge 3), and the global part of the "345" gauge 
symmetry. The fermion number "111" symmetry has an anomaly and the associated 't Hooft 
operator is, schematically, of the form (3_) 3 <9 + (4_) 4 (5 + ) 5 . 

2.2 GW kinetic terms and exact chiral symmetries 

To begin describing the proposal, we introduce three two-dimensional Dirac fermions, Vl/3, ^4, ^5, 
charged under the U (1) gauge group with charges 3, 4, 5, respectively. There is also a neutral Dirac 
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fermion, ^o- The fermion fields live on the sites, labeled by {x}, of a two dimensional lattice and 
their lattice action consists of kinetic terms: 



Skin — ^ ^ ^ g xD q x,y^q y j 
q=0,3,4,5 x,y 



(2.1) 



where D q is the Neuberger operator for a fermion of charge q, obeying the GW relation!! 

{D q , l5 } = D ql5 D q . (2.2) 

Here 75 is the appropriate two-dimensional matrix and the lattice spacing has been set, from now 
on, to unity. The lattice action (12, ip has a large number of exact global symmetries: 

[] (7(Vxf/(l) ?)+ , (2.3) 

9=0,3,4,5 

where C/(l) 9l ± acts only on the Dirac fermion field of charge q: 



(2.4) 



where P± = (l± 75 )/2 and P ± = (l±7 5 )/2, with 75 = (l-D)^. That 7 | = 1 follows from the GW 
relation f|2.2[) ; note that ^> q and ^ transform differently, which is perfectly natural in Euclidean 
space. The projector used for every \& 9 involves the appropriate Neuberger operator D q . 

That the symmetries in equation (|2.4p are all exact symmetries of the action (|2.ip follows from 
P^zD = DP±, yet another consequence of the GW relation (j2.2[) . The measure of integration, 
however, is not invariant under any individual J7(l) 9) + or U{l) q -. Instead, under an infinitesimal 
U{l) q) ± transformation (12. 4h with parameter a 9i ±, the measure changes as follows: 



n 

r=0,3,4,5 



n 

0,3, 

1 - ia qi± Tr (p± - P T 



U(l) q ,±: [[ d^> r d^> r 

r=0,3,4,5 



(2.5) 



Yl d$! r d$> 7 

r=0,3,4,5 



1 ±ia q , T -Tr D q j 5 



Eqn. (|2.5p implies that for vectorlike symmetries U(l) q v = a<7,-)> there is no Jacobian and 

thus they are true symmetries of the theory. On the other hand [5,31], since —^TxD q j^ = n+ — n°_, 
the difference between the number of left- and right-handed zero modes of D q , the continuum 
violation of charge for anomalous symmetries is reproduced by the nonzero Jacobian. 



2.3 Reduction of the global chiral symmetries by GW- Yukawa couplings 

To construct our candidate "345" chiral lattice theory, we introduce a unitary higgs field, 4> x , living 
on the lattice sites. We will use <p x to write all possible Dirac and Majorana Yukawa terms that 
violate all symmetries (|2.3p of the kinetic term (|2.1h except: 

?7(1)3,- x C/(l) 4) - x C/(l) 5 ,+ x C/(l) ,+ . (2.6) 

For a review of the GW relation and exact chiral symmetry on the lattice, see, for example, [8] and references 
therein. 
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To this end, we relate the Dirac fields ^ g to their chiral components: * g ,± = P±^q, *g,± = ^ q P^', 
note that the definition of the ^>± chiral modes is now both momentum and gauge-background 
dependent. We then write down the most general Dirac and Majorana couplings — expected to give 
mass of order the inverse cutoff — to the fields: 

X + = *3,+ *4,+) 

/ *5 



(2.7) 



where T denotes transposition (we treat unbarred Dirac spinors as columns and barred ones as 
rows) of the formH 

Smass = VY, X +* M ^x, <&) Y- x . (2.8) 

x 

Both Dirac and Majorana masses are to be included for the mirror fields: if Majorana masses are 
omitted, there will be extra unbroken chiral symmetries and unlifted zero modes in an instanton 
background, resulting in failure to reproduce the 't Hooft vertex. Therefore, Majorana masses are 
needed in order to introduce a violation of fermion number into the lattice, as has already been 
noted in [32]. 

Instead of writing explicitly the entire matrix M(cj) x , 4>*), which depends on powers of <j> and 
(p* determined by gauge invariance, as indicated in (|2.8p . we give an example of a Dirac mass 
term: ^ , o ,-(0*) 3| I / 3,+ + *3,+0 3l ^O,-> an d of a Majorana mass of the form: ^ , 5,-72</ )8 (^ , 3,+) T — 
^3 ' +J2(4>*) 8i &5,-- Here 72 is the hermitean gamma matrix that appears when Majorana masses are 
written using Dirac spinors, while 4> is the unitary Higgs field. Thus, the explicit form of M(<p x , cj>*) 
in (12. 8p contains appropriate powers of </> and 72-insertions. The general mass matrix (12.81) violates 
all [7(1) symmetries except the desired 1/(1)3 - x [7(1)4, _ x [7(1)5,+ x ^(l)o,+ symmetry from (|2,6p . 

We continue by noting that not all symmetries (|2.6f) of the lattice action (|2.8j) , p. Ill are respected 
by the lattice path integral measure. The measure is only invariant under three symmetries: the 
^(1)345 and the [7(1)133 chiral symmetries — linear combinations of [7(1) 3 ,_ x [7(l)4 i _ x [7(1) 5 ,+ 
with coefficients 345 and 133, respectively — and the [7(1)0,+, which acts only on the n+ = P+^o 
neutral field, whose dynamics is expected to decouple from the physics of the charged sector. The 
third linear combination of the first three [7(1) 's in equation (|2.6p — the fermion number symmetry 
of the light charged fields, which can be taken to be the "111" symmetry — has an anomaly exactly 
reproduced by the Jacobian, eqn. (|2.5p . of the corresponding transformation of the measure; see 
[5], [31], and references in [8]. Thus, the lattice theory obeys exact Ward identities, including the 
anomalous ones. For example, using (12. 5|) one finds that the 111 transform of an operator O obeys 
the exact lattice Ward identity: 

(5a lu O) = i\ (O TV [75(1)3 + D 4 - D 5 ))) . (2.9) 

The continuum limit expansion Ticy^Dg ~ J d 2 x e^F^ [31] implies that the anomalous Ward 
identity (|2.9p has the expected continuum limit. 



6 The presence of neutral mirrors is necessary to preserve two dimensional Lorentz invariance. 
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At the end of this section, it is worth noting that this proposal carries some of the flavoiu of an 
earlier construction of Eichten and Preskill [13], attempting to use strong four- Fermi interactions 
to decouple mirrors and doublers (it is clear that integrating out our short-ranged 4> x will produce 
strong multi-fermion interactions of the mirrors). Their proposal is known not to give rise to a 
chiral gauge theory (see [34], where the similarity with Yukawa models was also used). In our case, 
the modified lattice chiral symmetry that leads to exact decoupling of the chirality components 
only allows us to make use of the Yukawa analogue of the strong four-Fermi coupling symmetric 
phase (see the Appendix of [13]) — a phase with unbroken gauge symmetry, where all fermions that 
participate in the strong interactions are massive. 

2.4 Action, partition function, and dynamics 

To ensure that the dynamics of our lattice model reproduces that of the desired unbroken chiral 
gauge theory, we need to demonstrate the existence of a strong- Yukawa-coupling symmetric phase 
with chiral spectrum of massless fermions (recall again the strong coupling analysis of [22] which 
showed that in the waveguide model the spectrum in this phase was vectorlike). Remarkably, 
as we will find below, to leading order in the strong Yukawa coupling expansion and small gauge 
coupling — precisely the regime where the waveguide idea broke down — there appear no new massless 
modes and the spectrum of the unbroken gauge theory is chiral. 
The total action of the lattice model is, finally: 

S = S\yn son + Ski n + S mass + S K , (2.10) 

Skin is defined in (|2.ip , S ma ss — in (|2.8H , Swuson is the usual plaquette action for the link variables 
U x , x +fi, appropriately modified to restrict the gauge field path integral to admissible gauge field 
backgrounds, see [8], and S K is the action for the charge-1 unitary Higgs field: 

S - = \ E £ t 2 " ( & U **+t <^+A + h -c )] • (2.H) 

The dynamical issue that needs to be addressed is the existence of an "unbroken" phase where (j) 
is disordered (analogous to ((f)) = 0, versus ((f)) ^ 0, in four dimensions), such that the gauge boson 
is massless. 

In the case without fermions, it is well known [35] that theories with unitary Higgs fields 
(contrary to "everyday" continuum intuition) exhibit a symmetric phase, for small enough k. The 
essential ide^ is that for small k large fluctuations of the unitary Higgs field — or, in the equivalent 
unitary gauge, the pure-gauge fluctuations of the gauge field U — are not suppressed by the action 
(I2.11j) and hence their correlation length is of order the lattice spacing. Thus, integrating out 
the rapidly fluctuating Higgs fields results in renormalization of the gauge coupling plus a tower of 
higher-dimensional gauge invariant local operators which are irrelevant for the long-distance physics 
of the gauge theory. This is most easily seen upon integrating over the rapid fluctuations of (j) x , 
or equivalently, the pure-gauge part of U, by explicitly performing the strong-coupling (small k) 
expansion. The leading correction is a small, 0(k 4 ), shift to the inverse gauge coupling constant, 

7 We thank David B. Kaplan for pointing this out to us. We also note that a proposal to decouple the mirrors 
by combining the (approximate) lattice chirality of domain wall fermions with the Eichten-Preskill ideas was made 
earlier in [33]. 

8 Sometimes called the "FNN mechanism" [36]. 
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g 2 — > g 2 + k 4 (see section II. C. (a) of [35]), while the tower of higher-dimensional gauge-field 
dependent operators is suppressed by increasing powers of the small correlation length of (\> x - 

Since the gauge theory (|2,10p is asymptotically free, to study the continuum limit it is sufficient 
to begin at leading order in the g — ► expansion; here we will confine our analysis to this limit. 
The limit freezes the gauge degrees of freedom to U = 1; see also discussion in Section [331 The 
resulting theory is a unitary Higgs- Yukawa model, equivalent to the XY model coupled to fermions, 
whose phase structure can be studied in various limits. We are interested in the symmetric phase 
of the lattice XY model and will study k < k c while also taking y — > oo (cf. eqn. (|2.8jl ). 

In a trivial gauge background, the lattice partition function of the model (|2.10p factorizes into 
a product Z = Zu ght x Z mirror , with: 

Zlight = /n d *3,-^4,-d*5, + ^o,+ e- 5 «»^ i8 ' lt ] , (2.12) 

J X 

/~1 — r . _ - - - omirror ["iiv mirror] q \ jA c mirror ji 

J];^3, + ^4,+^5,-^0,-#e- S ^ I* ]-S K [4>]S mass l* ,*] j 
x 

where we used the decomposition of the Dirac fermion measure into left- and right-chirality mea- 
sures, d^f = dfy+d^-, see following Sections, and, for conciseness, omitted the conjugate fields 
in the measure. We also denoted collectively by ^u g ht the fields ^3,-, ^4,-, ^5,+, ^0,+, and by 
^mirror the heavy charged mirrors ^3,+, ^4,+, ^5-, and the neutral ^0— • The "mass" term is 
given by equation (|2.8p and the kinetic term for <f> by (|2.1ip . 

The most important point is the splitting of the kinetic terms (]2.ip into light and mirror modes 
in (|2.12p . This follows from the identity, which also holds in an arbitrary gauge background: 

* q D q * q = V qj+ D q -$ qj+ + #,,_£>,*,,_ , (2.13) 

where the cross-terms vanish due to the GW relation (|2.2|) . Thus the mirror and light partition 
functions factorize at g = 0. We stress that the GW relation was crucial in order for (|2.13p to hold; 
there is no other way to achieve (|2.13p . and hence the factorization (|2.12p . on the lattice, for D q 
free of doublers. 

Finally, let us study Z m i rror and its effect on the light modes, in the y — > 00 and k < k c 
limit. Of particular concern is the possible appearance of extra massless states and the associated 
vanishing of the mirror determinant. To this end, we redefine the mirror fermion fields in (|2.12p . 

^4,+j ^5,-1 and the singlet ^o,- 5 by 1/^/y- This multiplies their kinetic terms by 1/y. Thus, 
as y — > 00, the mirror fields kinetic terms vanish, and the mirror fermion action consists solely of 
a Yukawa term given by (12. 8h with y = 1: 

Smass = Y, X +- M (^' <&) Y ~ - ■ ( 2 - 14 ) 

x 

We can now perform the integral over the mirror fermions in Z m i rror . The notation in the above 
equation might lead one to believe that the mirror fermion path integral is a product of strictly 
local factors: 

J]detM(<^,0*) , 

X 

where M(4> x , 0*) is the 8x8 matrix from (|2.8l2.14p . which is gauge covariant and depends only on 
the values of <j) x at x. Since <f> x is unitary, this would imply that the gauge invariant detM((p x , </>*) 
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is a.) (^-independent and b.) nonzero. Hence, we would have been led to believe that the dynamics 
governing the fluctuations of the pure- gauge degrees of freedom is unaffected by the mirror fermions, 
to leading order in 1/y. Moreover, this argument would also imply that there is no fine-tuning, at 
large y, required in order to keep the XY model in its high-temperature phase. Finally, a constant 
determinant, as would be obtained at large-y from the above argument, indicates that there are no 
massless fermion states, as a massless fermion state is expected to lead to a zero determinant. 

The true story, however, is more complicated than the discussion of the previous paragraph. 
This is due to the fact that the various ^>± chiral components which enter X+ X ,Y- X (12. 7p and 
Smass are somewhat smeared due to the nonlocality of the chiral projectors that define the chiral 
components for the barred fields. However, the extent of the nonlocality of the ^-^-component 
is small, governed by the range of nonlocality of Neuberger's operator, which is of order of the 
lattice spacing with an exponential tail, as the analysis of [37, 38] shows. Hence, one expects that 
the qualitative arguments of the previous paragraph still hold, together with the conclusion that 
the mirror fermion fluctuations do not significantly affect the pure gauge fluctuations and their 
determinant is nonzero. Section [3.31 is devoted to verifying this conjecture. 



3 The simpler "toy" model 

3.1 Definition of the toy model: action and symmetries 

In our analytic and numerical study, we will use a simpler model that captures the main features of 
the mirror sector dynamics at g = 0. The model has a minimal field content, allowing an exhaustive 
study of the phase diagram using numerical methods with the computer resources available to us. 

Our toy model is a U(l) lattice gauge theory with one charged Dirac fermion, ip, of charge 1, 
and a neutral spectator, x- The desired spectrum of light fields in the target theory is the charged 
ip+ and the neutral x~- The chirality components for the charged and neutral fermions are defined, 
as in the previous section, by the projectors which include the appropriate Neuberger operators D\ 
or Dq for the barred components. The fermion part of the action of our toy model is: 

Smirror 

(3.15) 

Sught = (ip + ,D 1 i/) + ) + (x-,D X-) 

Smirror = At^-) + (x+, D 0X+) 

+ y{(^,f X +) + (x+,H>-)+h[(ip-,<t>l2X+) ~ (x+,720V-)]} • 

Here <p x = e irix is the unitary higgs field and we do not show its kinetic term as it is the same as in 
(|2.10p . The brackets indicate both summation over coordinates and an inner product of spinors, for 
example (^_, <j)*x+) = Y, i>- x<P x X+ x and Z?i^+) = Yl xDi xyip+y There are two Yukawa 

x x,y 

couplings in the model, y and yh, which are both taken real. The coupling h measures the ratio of 
the Majorana to Dirac mass, while y is the overall strength of the Yukawa coupling. The S m i rror 
term above is the analogue of (|2.8p in the "345" theory. 

When y = h = 0, the lattice action (|3.15p has four global U(l) symmetries, as every chiral 
component can be rotated independently, as in Section 2. When both y and h are nonzero, there 
are only two U(l) symmetries, acting on tp + and respectively. The first is the anomalous global 
part of the gauge group and the second is the global symmetry of the spectator fermion. When 
h = 0, the Majorana mass vanishes, and we have one extra exact U(l) that also acts on the charged 
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fields and leads to extra zero modes. Hence, we also need h 7^ 0, as we will explicitly see below; 
note that the phase diagram for a four-dimensional model very similar to ours, without Majorana 
mass terms has recently been considered in [28], using analytical methods!! 

The continuum target theory that f)3. 15H is conjectured to flow to is the chiral Schwinger model, 
an anomalous U(l) gauge theory with a single right-handed massless fermion of unit charge 
a model solved in [30,39]. Since our paper is confined to the study of only the mirror fermion-Higgs 
dynamics, without gauge fields and without taking massless fermion loops into account, we do not 
have to cancel anomalies of the light spectrum; minimality of the field content is why the "toy" 
model was chosen as an arena to study the strong mirror-sector dynamics. 

The fermion action (13.15p . the gauge action, and S K (the action of cf) in (|2.1ip ). combined with 
the usual Dirac fermion measure and the measures for the U(l) link fields and <j> x , completely define 
the lattice partition function of the toy model. As a matter of principle, the lattice model defined 
above could be, right away, subjected to numerical simulations including both dynamical gauge 
fields and fermions. However, this is beyond the scope of this paper, which is restricted to the 
g = case. Clearly, this is only a first relatively simple step in the full analysis of the dynamics 
and the utility of the proposal. Working out the dynamics at g = is also a litmus test for the 
proposal: if the analysis fails to yield a chiral spectrum or an unbroken gauge symmetry (more 
precisely — its global counterpart) at this point, the proposal has to be abandoned and no further 
study at g / is warranted. 

3.2 The toy model partition function at g = 

When g — > 0, only zero-action fluctuations of the gauge field contribute to the path integral. In 
other words, U xy = with u x = e lotx . The uj x fields appear in the kinetic terms of (|3.15p and 

in the Yukawa couplings, via the projectors used to define chirality components ip± of the charged 
?/>. We can remove the oj x dependence in the lagrangian by a local field redefinition of "if), if), and 4>: 

iP x = e* a * V4, 4 = & e~ ia * , <f> x = e-^^ , (3.16) 

with unit Jacobian. Thus the path integral over the tu x decouples, yielding the volume of the group. 
We are left, dropping the primes, with a nontrivial integral over ip, x, and (p. 

Now, as explained in Section 2, the partition function in a trivial gauge background splits into 
a product of mirror and light partition functions. As also explained there, for the purpose of the 
study of the g = dynamics of the mirror fermion sector, the chiral split of the measure should be 
worked out explicitly, in order to make the decoupling of the light and mirror sector manifest and 
to facilitate a Monte Carlo study of the mirror sector dynamics. 

To work out the ±-chirality split, it is useful to define the measure in terms of eigenvectors 
of the Neuberger operator. Since the chiral projectors involve the Neuberger operator, using its 
eigenvectors to define the measure facilitates an easy split of the measure into chiral components. 
The details are straightforward but somewhat tedious and are given in Appendix A. For the reader 
interested in the final formulae, the important expressions are eqn. (jA.30h — the expansion of the 
GW-chirality components of spinors in terms of GW momentum eigenstates — and eqn. (|A.17p . the 
final form of the measure relevant for its chiral split (also reproduced here in (|3.19j) . (|3.20p ). In 

9 We also need h ^ 1: at ft = 1 the mirror fermion spectrum has an exact zero mode at finite N for arbitrary <f> x 
backgrounds, see Section \3. 2 1 
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this section, we only present the final expression for the partition function and measure resulting 
from applying these formulae. 

We start with the partition function written in terms of position space variables, the charged 
{ip±x,i>±x} and neutral {x± x>X±a;}) in terms of which the chiral split of the measure is quite 
complicated. Then, we change variables, using (|A.30|) . in both the action and measure in terms 
of appropriately chosen — corresponding to the GW eigenvalues with positive imaginary part — 
momentum space variables {tp± x , ip± x } -> {a k ±,a k ±} and {x±x,X±x} -> {Ac±,Ac±}- After 
performing all the required summations over x, y, the results for the various components of Su g ht 
and Smirror are as described below. 

The kinetic terms for the mirror and light fermions from (|3.15p become, see fjA. 13[) : 

^mirror = E Ak ("k-"k- + Ac+Ac+) > ( 3 - 17 ) 
k 

S&u = E A k («k+«k+ + Ac-Ac-) ■ (3.18) 
k 

Here Ak is the eigenvalue of the free Neuberger operator with positive imaginary part, eqn. (|A.22|) . 
corresponding to momentum k. The most important properties of Ak are that Ak — 1 is a complex 
number of unit modulus, a consequence of the GW condition, and that Ak only has a zero at k = 0, 
see (IA.22p . The momentum sums here and below are over k±, = 1, N. 

The measure of the g = partition function, see eqn. (|A.17|) . splits naturally into a measure 
for the light fermions ak + , ak + , Ac - > Ac - : 

N 

Y[ #+ x dip+ x dx- x dx- x = ]7 16(1 - A k )(iak+(iak+(i/3k_d/3k- , (3.19) 

X fcl,fc 2 = l 

and a measure for the mirror fermions Ok <5k Ac +, Ac +'■ 

N 

Y\ dip- x dip- x dx+ x dx+ x = Yl !6(1 - A k )da k -(iak-cf/3k+d/3k+ • (3.20) 

X fcl,fc 2 =l 

To complete the definition of the mirror fermion partition function Z m i rror , we need to write 
the Yukawa couplings of the mirror fermions in S m i rror ()3. 15j) to the unitary Higgs field. To this 
end, we first define the Fourier transforms of (f> x = e ir,x , with ujn = e N : 

^4e ^ k ' x eivx > ^ kx e ' ivx = (*-*)• • ( 3 - 21 ) 

X X 

The two Yukawa couplings in the toy model lagrangian are, then: 

-SmtZr = ^E( 2 " A k) (ak-/3p + $ k _ p + 4 + «p-$k- P e^ p ^ k) 

y ^ , v 

k,p 

1 Q Ma i _ ( q 22^1 

yj^ mirror 
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where the momentum-dependent phase factors e lipk are defined in eqn. (|A.28p . 

The existence of an exact mirror- fermion zero mode at h = 1, alluded to in Section 13,11 and 
evident in all our numerical results, for arbitrary <j) x backgrounds, can now be shown. The mirror 
zero mode at h = 1 is a Majorana-Weyl fermion: this can be seen by letting all but the k = 
components of /3i c +,/3k+ vanish, taking /3o+ = — i/3o+j and showing that the mirror action (|3.22[) 
vanishes if h = 1. 

The equations presented in this section completely define the partition function. The decoupling 
of the light and mirror partition functions is manifest. The mirror partition function Z m i rror is 

defined by means of the fermion measure (|3.20p . the measure for the unitary Higgs field, Y\ J dr] x , 

x 

and the action 5 mirror = S^ rror + S^ rror + S m ^ ror + S K , with the various terms defined in (|3.17p . 
([32ZD, and (I2TTTD . 

As mentioned in the Introduction, the split of the measures allows us, from now on, to concen- 
trate on studying the properties of the mirror fermion-Higgs sector with partition function 'Zirnirr&r- 

3.3 Evidence for the existence of a symmetric phase at strong Yukawa coupling 

In this Section, we present our study of the symmetry realization in the mirror theory at strong 
Yukawa coupling. We remind the reader that the two possibilities are: 

1. A disordered ("symmetric") phase, where correlation lengths are of order the lattice spac- 
ing. The absence of order in the Higgs field is the analogue of the symmetric phase in four 
dimensions. The gauge boson would remain massless, apart from the mass generated by the 
Schwinger mechanism. 

2. A quasi-ordered phase, where correlations fall off according to a power law. In this phase the 
gauge boson would acquire a mass-squared of order Kg 2 . It is the two-dimensional analogue 
of spontaneous gauge symmetry breaking in four dimensions. 

In the quasi-ordered phase there are long wavelength modes that could mix with fermion bilinear 
composites and form light mirror sector states. Hence any indication of critical behavior in suscep- 
tibility data will be a sign of possible problems for the decoupling of the mirror sector. To this end, 
we perform a Monte Carlo simulation of the mirror theory, whose partition function was defined in 
Section [3. 21 all averages are computed using the Monte Carlo techniques described in Appendix ICl 
Apart from a few remarks below, all simulations are done in the limit when the kinetic term of 
the mirror fermions, eqn. (|3.17p . is neglected, i.e. in the strong Yukawa coupling limit y — > oo (cf. 
Appendix IC.ip . The ratio of Majorana to Dirac Yukawa couplings, the parameter h in (|3.22j) . is 
kept fixed. For N = 4, 8, 16, our results indicate that in this limit: 

• There is no sign problem for the fermion determinant at h > 1. We have no analytic proof for 
arbitrary 4> x background (we have been able to show positivity for a constant background and 
to first order in small fluctuations). However, based on our simulations of the determinant 
in many (tens of millions) arbitrary 4> backgrounds, not a single configuration has lead to a 
negative determinant for h > 1. The results of the simulations and the presence of the exact 
zero mode at h = 1 suggest that it should be possible to find an analytic proof of positivity. 
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For < h < 1 there is a sign problem and indications of critical behavior, as described below. 



• All correlation lengths, for k = 0.1, 0.5 (both smaller than the Xy-model k c ~ 1.12) are 0{a) 
and there is no critical behavior for h > 1 as discussed below. 

• Within errors, for all h > 0, there is no direction of r\ that is favored. The shift symmetry 
r] — > rj + const, seems to be preserved. However, this is not a test of quasi-order, as we need to 
look at the response to a small external perturbation — hence, we perform the susceptibility 
measurements. 

• The susceptibility, for the same values of n, shows a large response for h < 1 that increases 
with volume, indicative of quasi-order, and a small response for h > 1, indicative of the 
absence of quasi-ordering. 

• The vortex density decreases as h decreases towards h = 1, as in the Berezinski-Kosterlitz- 
Thouless transition towards the low-temperature "confined-vortex" phase. 

In this summary, and in the discussion that follows, it is important to keep in mind that the 
phase transition at h = 1 is not sharp at finite N. Therefore when we speak of properties at, say, 
h > 1, what is implied is behavior at h > 1 + e with e > a small number that tends to zero as 
N -> oo. 

3.3.1 Scalar field susceptibility 

We begin with a study of the scalar field susceptibility, which is a measure of the long distance 
correlations of the XY-model field <j) x = e irix . In the "high-temperature" disordered phase of 
the XY model the susceptibility is small, of order the lattice spacing or less, depending on the 
value of k, < k c . The critical coupling for the pure XY model on a square lattice is known to be 
k c ~ 1.12 [40]. In the "low-temperature" phase with algebraic long-range order, the susceptibility 
increases with the volume of the system. 

The susceptibility, which can also be thought of as the zero-momentum propagator of (fi x , is 
defined in the usual way: 



For an explanation of the meaning and numerical implementation of the (. . .) average, see Appendix 



We now present the results of numerical simulations of the scalar field susceptibility (|3.23[) . 
First, in Fig. [H we show the susceptibility for k = 0.1 for N = 4,8,16, as a function of h; the 
pure XY-model values are shown by horizontal dashed lines. We see that for h > 1 the fermions 
do not have an appreciable effect on the scalar field susceptibility. Similar results are obtained for 
k = 0.5, with the resulting overall rise in susceptibility, due to the closeness of k to k c = 1.12, and 
are shown on Fig. [2j Finally, in Fig. [3l we give the result for k = 0. It is qualitatively similar to 
the k = 0.1,0.5 results. 

The results of the simulations of the scalar susceptibility displayed here show that for h > 1 
the fermions have a negligible effect on the scalar field susceptibility. The results of this Section 




(3.23) 



C. 
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h 


N 


X 


N 


X 


N 


X 


0.01 


4 


3.15(7) 


8 


3.42(16) 


16 


3.42(14) 


0.005 


4 


3.15(6) 


8 


3.41(16) 


16 


3.41(14) 


0.001 


4 


3.14(7) 


8 


3.41(16) 


16 


3.40(14) 



Table 1: k = 0.1 results for small values h — > 0; note that no h — ► extrapolation is required, as the limit 
is already obtained within error and the susceptibility \ appears to approach a constant value as N — > oo. 



support the conjecture that at strong Yukawa coupling the theory remains in the symmetric phase. 
We also note that the susceptibility of the theory with fermions, for h > 1, is somewhat lower 
than the susceptibility of the pure YY-model with the same value of re, for iV> 1. This can be 
interpreted as a small renormalization of re to smaller values due to the fermions, which appear to 
push the theory further into the symmetric phase. 

On the other hand, for h < 1, there are indications of critical behavior and the scalar suscepti- 
bility grows with the volume of the system, as appropriate for the low-temperature phase of the XY 
model. Also, for h ~ 1, the fermion determinant vanishes and there is a massless mirror fermion 
state (see the footnote in Section 13.11 and the numerical analysis of the spectrum in Section I3.4|) . 
Clearly, at these values of h, the sensitivity to y < oo is enhanced, see Section lC.il in the Appendix 


Finally, for h — > 0, in Table [1] we show the re = 0.1 results for the Higgs susceptibility. Clearly 
there is a finite h — » 0, N — * oo limit, contrary to the finite-size scaling that would occur in a quasi- 
ordered phase. We conclude that the Higgs is in the disordered phase also for h = 0, in agreement 
with the analytic [28] and numerical [29] results in a similar model (h = 0) in four dimensions. 



3.3.2 Binder cumulant 

The Binder cumulant [24] is a quantity probing higher-order correlations, defined as: 

U = 2-<^%, (3.24) 
(|M| 2 ) 2 V ; 

where M is the total "magnetization:" M = 52 x <f> x , \ M ? = MM*, |M| 4 = (|M| 2 ) 2 . The 
definition (|3.24p is chosen such that in the pure YY-model U interpolates between and 1 as the 
temperature decreases from infinity to zero: 

lim U = 0, lim U = 1 . (3.25) 

This is easily seen by noting that for k = 0, there is no intersite correlation and hence U can be 
computed by treating the <p x as random, delta-correlated, field leading to U = 0. On the other 
hand, for k = oo, there is perfect order in each configuration and the <p x are frozen, giving rise to 
(7 = 1. 

We show the results of the measurement of the Binder cumulant for two values of k: re = 0.1 
in Fig. S] and re = 0.5 in Fig. [SJ In each case, as a function of h the Binder cumulant behaves as 
if h < 1 corresponds to the broken phase of the XY model and h > 1 — to the unbroken phase, 
consistent with the measurement of the scalar susceptibilities of the previous section, and thus 
providing further evidence for the persistence of the symmetric phase at strong Yukawa coupling 
and h > 1. 
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3.3.3 Vortex density 

The Berezinski-Kosterlitz-Thouless transition in the XY-model can be interpreted, see [23], as due 
to the deconfinement of vortices above the critical temperature (at k < k c ). The definition and 
algorithm of finding the vortex density — the number of vort ex- ant i vortex pairs per XY spin — is 
described in [41]. 

We studied the vortex density as a function of h. We show the results, for k = 0.5, on Fig. [6j 
It is clear from the figure that the vortex density of the deformed XY model at h > 1 is slightly 
higher than the value appropriate to the pure-XY model for n = 0.5; this is consistent with the 
susceptibility measurements, where the small effect of the fermions at large y is to push the theory 
further into the symmetric phase. 

To summarize, in Sections 13.3. H 13.3.21 and 13.3.31 we presented measurements of the pure scalar 
probes of symmetry breaking — the scalar susceptibility, the Binder cumulant, and the vortex den- 
sity. The results give strong support in favor of the existence of a symmetric phase at strong 
Yukawa coupling, at h > 1, K < 1. 

3.3.4 Fermion bilinear susceptibilities 

In this Section, we consider another probe of symmetry breaking — the fermion susceptibilities. As 
discussed in the Introduction, a scalar <p x and a charged mirror fermion, can produce a neutral 
fermion bound state at strong coupling. Similarly, one expects that a charged and a neutral fermion, 

and x+7 can bind into scalar bound states with the quantum numbers of 4> x . 

The two composite charged complex scalars of lowest dimension that we can construct out of 
bilinears of the mirror fermions are rp— x+ and ■0-72X+- We then define the corresponding suscepti- 
bilities; similar to the scalar case (13.230 they can be interpreted as the zero- momentum propagators 
of the corresponding composite scalar fields. Thus, we introduce the "Dirac" susceptibility: 

XF = ^2(tp- x X+x X+ y^-y) , (3-26) 

X 

and the "Majorana" susceptibility: 

Xf = ^(''P- xl2X+ x X+y72^- y ) • (3-27) 

X 

In four dimensions a disconnected part would also have to be subtracted; however, the disconnected 
part always vanishes in two dimensions due to the absence of spontaneous symmetry breaking. 
The explicit form of xf and x'p 1X1 terms of averages over the variables of integration introduced 
in Section [3.21 is given in eqns. (1B.16[ IB,19j) of Appendix [Bl Note that the dimensionless fermion 
susceptibilities scale as Y -2 , at large dimensionless Y = ya. 

The results for the modulus and argument of the Dirac susceptibilities are given, for k = 0.5 
in Figs. [7] and respectively. Similarly, for the Majorana susceptibility, the results for k = 0.5 
are shown in Figs. [12] and [T3l respectively. The same plots of the magnitude and arguments of the 
Dirac and Majorana fermion susceptibilities, this time for the smaller value of k = 0.1 are given in 
Figs. [H [101 QTJ Q3] and lead to the same qualitative conclusions. 

We observe that the behavior of their modulus follows, as a function of h, the behavior of 
the scalar susceptibility and the Binder cumulant. We note the following dependence on system 
volume: 
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• For k = 0.1 and h > 1, the fermion susceptibilities never show a trend of increasing with N. 
While the curves for k = 0.5 are not monotonic, there we also find no evidence for finite-size 
scaling that would indicate the presence of light modes. 

• x'f shows a trend of increasing for h < 1, for either k. For xf, the trend depends on k. 
This indicates that the two susceptibilities are indeed independent probes of long distance 
behavior. The rise in either with N is indicative of light modes that exist at h < 1. 

The numerical results for the phases of the Dirac and Majorana susceptibilities at h > 1 are 
consistent with phases equal to and ir, respectively. 

As a matter of comparison, we have also studied the fermion susceptibilities in the broken 
phase, upon increasing k > k c . In the broken phase, the fermion susceptibilities (as well as the 
scalar susceptibility) show a dramatic increase with the volume of the system, which is clearly 
absent at k < 1. 

Overall, our results for the fermion susceptibilities show no sign of critical behavior at h > 1, and 
thus provide more evidence for the existence of a strong Yukawa symmetric phase. 

3.4 Evidence for the absence of massless mirror fermions at strong Yukawa 
coupling 

In this Section, we turn to the fermion spectrum at strong Yukawa coupling !/> 1 and k < 1. We 
study two kinds of correlators: those that probe the neutral fermion spectrum and those probing 
the charged fermion sector. All correlators that violate the U(l) symmetry, i.e. of the form t/> + x_, 
vanish, due to the absence of spontaneous symmetry breaking of a continuous symmetry in two 
dimensional [19]. 

In the disordered phase the mirror fermions obtain mass by the mechanism explained in the 
Introduction [18], [13-17]: a charged fermion can bind with the charged scalar and pair up with 
the neutral fermion into a heavy neutral bound state. It is also clear that in the disordered phase, 
the charged mirror fermion spectrum (if charged stable single-particle states exist) is vector-like: 
the neutral mirror fermion can bind with the scalar into a charged right-handed fermion, which 
then can pair up with the charged left-handed fermion into a heavy charged state. We will study 
correlators that probe the spectrum of both types. 

Studying the fermion propagators in the y — > oo limit will be a reliable guide to the spectrum of 
the mirror fermion theory at y 3> 1, so long as there are no massless states that invalidate the 1/y- 
expansion and make the y < oo corrections important. Further comments on the y < oo limit are 
given in Appendix lC.lt especially important is the numerical result InPf A4 = IndetM mentioned 
there. Thus, the absence of massless states for h > 1 that we find ensures the self-consistency of 
the analysis. To probe the fermion spectrum at large y, we thus study the Green functions which 
do not vanish in the large- y limit (note that at infinite y, all +- (- and correlators vanish). 

To study the fermion spectrum of our model, we begin by first considering the charged fermions. 
Using ?/>-, X+; an d <fi, the three simplest unit-charge local fermion operators can be defined: 

Vix = ^-x , V2x = X+ x<t>* x , Vsx = X+ x <t>t ■ (3-28) 

10 As a check on our numerical simulations, we have verified this property in both the disordered and quasi-ordered 
phases. 
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Negative charge fermions can, similarly, be created by the local operators: 

f itx = i ; -x, ffe x = X+ x<t>x , V3x = X+x4>x- (3.29) 

The choice of this operator basis is motivated by simplicity and by the fact that the corresponding 
propagator (|3.3U|) can also be easily computed in the quasi-ordered phase, see eqn. ()B.7|) . and the 
change of the spectrum observed clS EL function of k. 

The propagator of the charged fermions is the connected two-point function: 

D?-(k) = J2 w N HX ~ y) ( Vlx f) c jy ) , i,j = 1,2,3. (3.30) 

X 

Diagonalizing Z??-(k) by a (bi-)unitary transformation gives rise to a diagonal matrix with entries, 
which we denote by l///j(k)(i = 1,2,3), where the //j are the eigenvalues of the inverse charged 
fermion propagator. A massless mode would yield a zero of some of the eigenvalues \ii for some k. 

The study of the spectrum is simplified at infinite y, when all ++ and Green functions of 

the mirror theory vanish (there is an extra symmetry of the partition function at y — ► oo forbidding 
such nonzero expectation values). Then it suffices to consider a reduced propagator, instead of 
(|3.30p . which we call S c , which we obtain by keeping only the nonvanishing entries of (|3.30[) : 

qe ( (^-yX+x^x) (X+ y<t>*y x) \ /oo-iN 

Clearly, the propagator S c obeys TrD c ^D c = TrS c ^S c in the infinite-y limit. Our Monte-Carlo 
simulation calculates the quantity f2( 2 )(k), defined through the last equality below: 

where the Fourier transform is defined, similar to (|3.3U|) . as S c (k) = Yjx^n*^ '^y-x' Eqn. (|3.32p 
shows that a massless particle will manifest itself as a zero of Q( 2 ) at some k. To look for light 
mirror modes, we will simply plot the value of minkfi( 2 )(k) as a function of h and the lattice size 
N, for y = oo. 

Next, we consider the operators creating neutral mirror fermions: 

Vix = X+ x , V2x = X+ x , V3x = M- x , rflx = x , (3-33) 
with their conjugates defined as in (13.291) : 

Vlx = X+ x , »)2 n x = X+ I > V3x = <t>* x i>- x , f]tx = <hrf>- x > (3-34) 
and neutral fermion propagator: 

Dg(k) = J2 U N HX ~ V) < Vix V? v ) , hJ = 1,2,3,4. (3.35) 

X 

In the y — > oo limit it suffices to consider the reduced neutral propagator, S n , keeping only the 
nonvanishing entries of f)3.35j) and combining them into a matrix similar to (|3.31l) : 



{X+y<PxV- x ) {X+y4>*x^- 



y - X ~\ (Xly^lx) {Xly^-x) ' ' 
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obeying TxD n ^D n = 2TrS n ^ S n in the infinite-y limit. We define the quantity £2W(k) by an equation 
identical to (I3.32p . but with S c replaced with S n ; our plots seeking to establish the presence or 
absence of massless neutral modes will similarly show the value of minkf2( 2 )(k) as a function of h 
and the lattice size N, for y = oo. Explicit expressions giving S c and S n in terms of the integration 
variables of Section \3. 21 are given in (|B.2l lB~7T2"j) of Appendix [Bj We use the expressions given there 
to calculate the O^'' 2 ) in our Monte-Carlo simulation. 

In the broken phase, when <j) x = 1 to leading order in perturbation theory, the form of S c (and 
S n , which coincides with it in this limit) and the perturbative spectrum, valid when k — > oo, are 
worked out in Appendix[B[ An important check on our Monte Carlo simulation is that, for k>k c 
it reproduces the perturbative spectrum of eqns. (|B.9l IB. lOh . including all the small eigenvalues 
found analytically to leading order in perturbation theory, see eqn. (IB.lOp . Furthermore, the small 
deviations of the Monte Carlo results from the perturbative spectrum can be seen to scale as the 
expected perturbative corrections to the tree-level result arising from spin-wave exchange, ~ /i 2 k _1 
(we find that these corrections become numerically significant for h > 2 for the smallest eigenvalues 
of a(k) _1 , see (IB^jl ). 

The quantities nrin^/H and mink^S provide a lower bound on the smallest eigenvalue of the 
neutral and charged inverse propagators, respectively. Finally, recall the large-?/ scaling (and that 
y has dimensions of mass) $7k ~ V- Introducing a dimensionless Yukawa coupling Y = ya, the mass 
scales as Ya" 1 . The continuum limit is N — ► oo, a = L p h ys /N, with fixed system volume L h . 

In Figure [T5l we display iVmink — the lower bound on the minimal value of the neutral 
fermion eigenvalue — as a function of h for k = 0.5 and N = 4, 8, 16. Figure [16] displays the 
corresponding quantity for the charged fermions. In Figures [T71 and [TBI respectively, we show the 
lower bound on the neutral and charged eigenvalues for k = 0.1. In units of L~^ ys , the lower bound 
on the mass can be obtained by multiplying the plotted value of by the dimensionless large Yukawa 
coupling Y. 

The numerical results for the lower bounds on the neutral and charged fermion propagator 
eigenvalues, /Uj, displayed in these figures, imply that: 

• The lower bounds on both the neutral and charged eigenvalues fa, for large Y, are ~ 
YN a L~^ ys , with a k- and (more weakly) /i-dependent exponent a > 0. This presents 
evidence that the mirror fermions decouple from the infrared physics in the TV — ► oo, fixed 
Lphys limit. 

• It is clear, by comparing Figs. [T5l and 1171 that the exponent a is an increasing function n. To 

1 2 

better study this dependence, we have also measured O k ' , for h = 2 and k = 0.6, 0.7, . . . , 1.4, 
in addition to values that we display. We found that in the interval 0.6 < k < 0.8 the 
exponent a ~ 1 (for larger values of k appropaching k c > 1 the susceptibilities rise with iV 
and the spectrum begins to approach that of the broken phase, see comments below). Thus 
for k in this interval, the mirror fermion masses are cutoff scale times Y. 

• If the lower bound on the fermion eigenvalues is actually saturated, one can conclude from 
the data that the charged fermions are generally heavier than the neutral ones. 

In summary, our results present evidence that there are no light states in either the neutral or the 
charged fermion mirror sector. 
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Finally, to compare with the fermion spectrum in the broken phase, the dashed lines on each 
figure show the minimum values of f2 for N = 4, 8, 16, calculated to leading order in perturbation 
theory (strictly valid when k — * oo) in the broken phase. This was done by using the broken 
phase propagator of eqn. (|B.9|) . As already mentioned, these perturbative results are in excellent 
quantitative agreement with Monte Carlo simulations in the quasi-ordered phase (performed for 
k = 10 3> n c — 1). The dips in the dashed lines correspond to the approximate solutions of (IB.lOj) . 
leading at finite N to almost massless modes Our results, discussed above, show that in the 
symmetric phase the corresponding massless modes are absent. 
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Appendices 

A Expansion of spinors in terms of chiral GW eigenvectors and 
chiral split of the measure 

In this Appendix, we derive an expansion of spinors in terms of Neuberger's operator eigenstates. 
This is useful in writing the chiral split of the measure in Section 13.21 We denote the Neuberger 
operator by D, obeying {75, D} = Dj§D and 75-D75 = D>. Recall that the eigenvalues A of D obey 
A + A* = AA* and thus lie on a unit circle in the complex plane centered at 1. For every eigenvalue 
with Im A / there is another eigenvalue A* whose eigenvector, a 2-component spinor, is = 
75*a- Here = X^\. The inner product notation we use is (\l/t^, ^\) = Yl x x = ^a',a- 

Orthogonality of eigenvectors with different eigenvalues follows from hermiticity of 75 D. Chiral 
components for the unbarred fields are defined in the usual way, *$>± = P±^f, while for the barred 
*, ^± = , where P ± = (1 ± %)/2, with 75 = (1 - D)~/ 5 . 

To define the Dirac fermion measure, we first expand an arbitrary spinor configuration in terms 
of the eigenvectors ^\ x of D as follows l 12 l 



1 We caution against concluding from the plots that the broken phase massless modes indicated in (|B.10[l are not 
present in the continuum limit for some values of h — it is easy to check that as N becomes large, the dips completely 
overcome the JV enhancement from the scaling of y, as the eigenvalues become closely spaced on the unit circle. 

12 Here and below, there is an implicit sum over multiplicities associated with a given eigenvalue A. 




(A.l) 



ImX>0 
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where c\ are Grassmann numbers^ The measure of the ^ integration is, therefore: 

Y[dV x = Yl dc x dc x , . (A.2) 

x lm\>0 

The contributions of the real eigenvalues A = 0, 2 are to be added both to (|A"1|) and (TO) . We 
note that there is no Jacobian in eqn. (|A,2j) . This follows from the fact that the GW operator is 
normal ([-D^D] = 0) and hence can be diagonalized by a unitary transformation; the Jacobians 
from ^ and ^> transforms cancel, for any gauge field background. 

Our next goal is to rewrite the measure (|A.2|) in terms of chiral components. This is easily 
accomplished for the unbarred field ^> x (keeping in mind that = 75 ^f\). Explicitly, we rewrite 
(fA~Tj) as follows: 

~C X + C A * *A*^ CA - C A * ^A x ~ ^ X* x ' 



*x 



E 

ImX>0 



V2 



V2 



V2 



V2 



ImX>0 



where we defined: 



a X ± 



A± x 



CA ± CA* 

^Ax±^A*x 
V2 



(A.3) 

(A.4) 
(A.5) 



Now we can use the relation (|A.4p to express the measure (|A.2p in terms of the variables a±, with 
a Jacobian that is only a numerical factor: 



17 £^3= J [ dc\dc\* = Y\ ^dax-da x + 



(A.6) 



7mA>0 



ImX>0 



and the expansion (|A.3P to express the action in terms of the chiral integration variables a± . The 
zero and cutoff (A = 2) modes are easily added to (IA.3j) . (|A.6h . as shown below. 

Now consider the more interesting ^ case. Note the following chain of relations: D^f\ = A^a 



* A A* 



* A A* 



A*^ A »7s = X*^\, where the 75 



hermiticity of D was used. We continue by expanding the row spinor ^ in terms of a complete set 



of \1/ A solutions as follows: 



(A.7) 



JmA>0 

again omitting the contribution of the A = 0, 2 eigenvalues that should be always added in ones 
mind. Similar to (|A.2h we define the barred fermion measure: 



Yl dv x = Yl dc ^ di 



C A * 



(A. 



ImX>0 



Here, we have displayed the contribution of only the eigenvalues A / 0,2; the contribution of the eigenvectors 
of the real eigenvalues A = and A = 2 are implicitly assumed present and will be added later. The reason is that 
they are always (anti)chiral in the usual sense: for A = 0, 2 the GW chirality is equal, up to a sign, to the usual 75 
chirality. In the free case of zero gauge background we will explicitly find these eigenvectors and see that the A = 
ones (one + and one — ) corresponds to k — {N, N), while the A = 2 ones (one pair of ± each) are obtained for the 
three values of momentum k = {N/2, N), {N, N/2), {N/2, N/2). 
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The set of relations above (|A.7|) are useful because they allow us to rewrite (|A.7[) in terms of P± 
components as follows: 



lm\>0 



*U=F*I. S (1-A) , , ^ k t*L(1-A*) 



CA 7T^ 1 + CA 



2 " 2 



v c A =f (1 - A*)c A . ^L=F^ a (l-A) 



lm\>0 

where the first line was obtained by applying the projector to (|A.7p and in the second we used 
|1 — A | = 1. Thus, similar to the definition of ct± of (IA.4j) . we can define: 

a A ± = —/= (c A =F (1 - A*)c A *) 

*a±- = * A . a (l -A)) , (A.10) 

We see from the above equation that the Jacobian of the transformation from c to a± variables 
depends on A; explicitly, eqn. (|A.8|) becomes: 

l[d^ x = J] dc A dc A *= J! 4(1 - A*)da A _da A+ . (A.ll) 

i lm\>0 lm\>0 

As a simple demonstration of the self-consistency of this procedure, we check that in the vec- 
torlike case we can use the measure in the final form (|A.6h and (jA.llh to calculate the determinant 
of D. To this end, write first the expansion of in terms of a±: 

q A+ + q A _ f «A + ~ «A_ t fA12 N 

lm\>0 v v v ; 

(recalling that zero and cutoff modes should be added as well) also note that the A* — 1 factor in 
the denominator does not vanish since A = 1 is not in the spectrum. Then, substitute (|A. 12|) as 
well as the a± expansion of Vl/ of eqn. (|A.3P into the GW action to obtain, using the orthogonality 
of the &\ wavefunctions and recalling that A = A*/(A* — 1): 

(#,£>¥) = ^* X (D*) X = Yl aA+aA+A + a A -a A -A. (A.13) 

x lm\>0 

Now recall the a± measure (jA.lip and the similar expression (|A.6|) (with constant Jacobian) for ^ 
to find that the path integral over the A ^ 0, 2 eigenvalues reduces to 

[l[d^ x d%e^ D ^ = c J] A 2 (l-A*) = c' J] AA* 

x lm\>0 lm\>0 

= c' J] A = c'det( J D)| V02 , (A.14) 

A^0,2 

i.e. the determinant of D, as appropriate (the factor d hides various numerical constants as well 
as the contributions of the A = 0, 2 eigenvalues). The point of this exercise was to show that the 
Jacobian factor in (jA.llj) was crucial to obtaining the right result. 
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We end this section by summarizing the important formulae that were derived: the expansions 
of the fermionic fields in terms of generalized chiral and antichiral Grassmann amplitudes a± , a± : 



E 

lm\>0 

E 

lm\>0 



V2 



V2 



a x+ + a X - , »A+ ~ ol\- 

Xx V2{\*-1) X * x . 



V2 



(A.15) 



(A.16) 



and an expression for the measure, combining (|A.6[) and (|A.11|) . in terms of these Grassmann 
amplitudes: 

Y[ d^f x d^ x = c ] | da\ + dax-dax+dax-(l - A*) . ( A -17) 

x lm\>0 

This defines a split of the measure in terms of chiral and antichiral fields. The zero and cutoff 
modes (the real eigenvalues of D) should be added to (|A,15l lA.161 IA,17|) . 

The usefulness of (|A.15l IA.16|) is that if we write the lagrangian in terms of ^-t, ^>± compo- 
nents, e.g., as in (|3.15j) . we have a straightforward way to express the corresponding terms via the 
integration variables ay-. For example, if a term contains x , we simply need use (|A.16|) and put 
all a_ to zero in that expansion, etc. 

A.l Explicit form for the GW eigenvectors for g = 0. 

Here we work out the explicit form of the eigenvectors of the Neuberger-Dirac operator in the free 
case. We begin by writing the free Neuberger-Dirac operator in momentum space: 



N 



u% s $>{k) , lo n = e 2 N , 



N 



(A.18) 



fcl,fc2=l 



N 



(A.19) 



fci,fe=l 



We also define = s'm^k^/N) and = cosfak^/N) and the following functions of kn = 1, . . . , N, 



a k = 
6k = 
Ck = 



1 - 2s\ - 2s| 

2s 2 C2 



V 7 ! + &44 
V 7 ! + 844 



(A.20) 



in terms of which D(k) takes the explicit form (we use 71 = (71,72 = 02,75 = 03): 

m= ( a k fc k + 6 k 
V «c k - 6 k a k 



(A.21) 
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The GW condition is equivalent to the relation a£ + 6^ + c k = 2ak, which is easily verified from 
(IA,20h . The eigenvalues of -D(k) are: 



A k = a k ± iy/bl + 4 . (A.22) 

With the help of this equation it is easily seen that the only real eigenvalues occur for k = 
(N,N), (N/2,N), (N,N/2) and (N/2, N/2). The (N,N) eigenvalue is the zero mode A = 0, while 
the other three correspond to the cutoff value A = 2. For every one of these eigenvalues there are 

two eigenvectors of opposite 75 chirality; these can be taken simply to be represented by 



1 



and ( J ) , as will be done explicitly in eqn. (|A.27|) . 



The eigenvectors with ImA > are explicitly given by: 

1 t fi* ( J 
V2N N \ 7^ 



^A k x = -i— uj k / ( *b k +c k ) , (A.23) 



*V*V) S E*L.*V,=^. (A. 24) 



and are easily seen to obey: 



a property assumed in deriving, e.g. (|A,13t [A.14j) . In particular, the expansions (|A.151 lA.161 IA.17]) 
hold by simply replacing the there by (|A.23|) : note that the sums and products there are over 
all k apart from the four points corresponding to real eigenvalues. 

Now consider the A = 2 eigenvalues; the A = ones can be considered in exact parallel. We 
proceed simply by recalling that A = 2 eigenvalues are always (anti)chiral in the usual sense. More 
precisely, for the unbarred spinors, GW and usual chirality coincide for any A, while for the barred 
spinors, GW and usual chirality are opposite for A = 2; this simply follows from the fact that for 
A = 2 we have 75 = (1 — D)j5 = —75. In our basis of gamma matrices definite 75 chirality &± 
spinors have the form: 

•♦-(;)••-(!)• (A.25) 

while definite 75 = -75 chirality * ± are V± = ^fP T = ^(1 =F 75)/2| A=2 = *(1 ± 7s)/2. Thus, the 
ty± spinors with A = 2 have the form: 

#+ = ( c ) , L = (0d). (A.26) 

This allows us to write the A = 2 solutions for two spinors and \& + (we use the fact that A = 2 
solutions in the free case occur for k = (N/2, N/2), (N, N/2), (N/2, N) only): 

, = 1 + J^ +Nx ^_ + u^+W^V) ( J ) (A.27) 



1 ( (N/2)x 1+ (N/2)x 2 (N/2)x 1+ Nx 2r Nx 1 + (N/2)x 2 \f 1 
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where the a, (3, 7 are the corresponding Grassmann integration variables of the A = 2 modes; also 
clearly the various position factors are either unity or of the form (— l) Xl , etc., which will be 
important when performing the summations after the further insertion in the action. 

Now, we put all formulae together, using the momentum space eigenvectors of -D(k) found 
above, to obtain the final form of the expansion of 'J, ^ in terms of GW eigenvectors. To this end, 
we first define the phase factor from the eigenvectors (|A.23h . now extended to include all k: 

f -%±£*= if k^(N,N),(% ,Af),(iV,f),(f ,§) 
e «Vk = J ^ v ' M ' / (A.28) 

\ 1 if k = (iV,iV),(f,iV),(Ar,f),(f,f) 
and introduce the unitary matrices and Vk: 

Then, the final result for the momentum-space expansion of the GW-chirality components of the 
Dirac spinors is: 

^ N / n 

k=i v 

■\ N ( 1 \ 

- £ «k + uft* (J ) , (A.30) 
k=i v 



1 - 

k=l 

1 - 

*+* = n E"k+^ kx (0 l)V k . 

k=l 

We note again that the barred spinors are treated as two-component rows, as explicitly evident 
in the above expansions. The sum in (|A.30P is over all momenta and thus includes also the real 
eigenvalues of the GW operator. Note also that Uf~ = cr^V^e^^a^ and that the chiral projectors 
for the *f? fields can be written in terms of U^, V^: 

P±k = ul P± U k = V£ P± V k (A.31) 

In this notation, the fact that as k — > the usual and GW chirality coincide is evident — in this 
limit, the unitary matrices U k , become the identity. 



B Fermion Green functions 
B.l Neutral fermion Green functions 

As explained in the main text, see eqn. (|3.36p . we define the following fermion observable, appro- 
priate in the large y limit, when only correlators between left and right fields are nonzero: 

on _ ( (X+ x <Py1p- y) (X+ x y) \ , R -, s 

X - y ~{(X T + xM T - y ) (Xlx^-y) ) 
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These correlators probe the spectrum of neutral fermions — as explained in the Introduction, the 
charged fermions ip- can bind with the charged scalars into neutral fermions, which then pair up 
with the neutral mirrors to form massive states; see the strong coupling expansions in [13, 14] for 
discussion of analogous phenomena. 

Then define, for every value of k, the four-by-four matrix propagator in momentum space, and 
use the expressions (IA.30P (x+ is obtained by replacing a + — ► /? + in the corresponding equation in 

dS35D ): 



E 



<A-kOp-k*-p> V^P.e^ (^_ k S p $;_ k ) V T k a_U p 



Here, ()A.29|) defines the U and V matrices and e llp is defined in (|A.28|) . We will study numerically 
the following quantity: 

a£> = 1 , (b.3) 

which, as discussed in Section 13,41 provides a lower bound on the smallest neutral fermion inverse 
propagator eigenvalue. 

It is instructive to consider the form of the neutral Green functions for constant <j) x = 1 back- 
ground, where they can be evaluated explicitly. This is, of course, only relevant for large values 
of k where the fluctuations of 4> x are nearly frozen. Note that for <p x = 1 there is no difference 
between the charged and neutral Green functions (defined in following section). The propagator 
Sk then takes the form: 

n_f (/3k«-k)^- k (f3 k ^)P + U k \ 

The biunitary transform of s^.: 

*-(!*)*(;«)■ 

takes the even simpler form: 

and makes explicit the fact that the rank of the matrix is two, the number of propagating complex 
degrees of freedom. The relevant expectation values entering the Green functions can be easily 
computed for constant (j> x and are given by (it is important to take into account the measure factor 
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1 — A k (|3.20|) when calculating the expectation values): 

</? k «- k > = - lh e-^ A ~ A \ + 2Xl ri(k)- 1 (B.7) 
2 -A? 



-i 



where 



</9-kOk) = -i/«^ k (l-A k )d(k) 
</3_ k a_ k > = -^-^ d(k)- 1 , 



d(k) = sin 2 % + /i 2 cos k , A k = 1 + e iek . (B.J 



2 

Combining (|B,6l IB.7|) we find: 

rf^f-^^^ ^+ ^ (B9) 

d(k) ^ -£^k e ^- k P_ -ih e^ k (l - A£) a- J 

The rank of s k is two; clearly, the nontrivial part is obtained by simply replacing a± and P± 
by unity. After the elimination of two rows and columns of zeros, the determinant of s k equals 
e^-^dQa) . Clearly, the denominator (jB.8[) of admits small eigenvalues for values of k where 
(assuming h > 1/2): 

| sin ^| = ^^ = . (B.10) 
1 2 1 V2/J 2 - 1 

For generic values of /i this equation has no exact solutions at finite N, except for h = 0, where the 
$k = corresponds to the three unlifted modes with A k = 2, and h = 1, corresponding to # k = 7r, 
where the k = (N, N) mode is unlifted. Other values of k for which (|B.10|) approximately holds 
occur, for sufficiently large h near k ~ 5, and give rise to small mass eigenvalues in the broken 
phase. 

As explained in the main text below (13.221) . the h = 1 exact zero mode persists for small k as 
well. It gives rise to the critical behavior seen in our numerical analysis at h ~ 1, and discussed in 
the main body of the paper. 

As a consistency check on our numerical simulations, we have also verified that for large values 
of k, in the algebraically ordered phase, where perturbation theory leading to (|B.9|) is valid, the 
mass matrix, including the small eigenvalues due to (|B.10|) . are reproduced by our Monte-Carlo 
simulations. 

Most importantly, however, as described in the main text, in the disordered phase with rapid 

1 2 

fluctuations of <j> x , our results for f2 k ' show that there are no small eigenvalues in the h > 1 regime. 
B.2 Charged fermion Green functions 

The correlators that probe the charged fermion spectrum, eqn. (|3.3ip . where charged bound states 
of the neutral mirrors x+ an d the scalars can pair up with the charged mirrors to form heavy 
fermions, appropriate to the large-?/ limit are as follows: 

Sy-*=( n' XXlv i V \ { ?-T X 5f y \) (B.H) 
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For every value of k, the four-by-four matrix propagator in momentum space is, then, similar to 
4 = »v k "- y (B.12) 

X 

(a k /3 p $_ k _ p ) e'«<7_ (/3pak* k _ p > F+C k 



p 



The quantity that we numerically study is defined as for the neutral propag ator (lR3l) : 

0® = , (B.13) 

V Tr v s k 

and similarly provides a lower bound on the smallest eigenvalue of the inverse charged-fermion 
propagator. 

B.3 Fermion-bilinear susceptibilities 

There are two composite charged complex scalars that we can construct out of bilinears of the 
mirror fermions — x \+ x and tp'E. x l2X+ x- Consider the corresponding susceptibilities: 

&x-y = {4>- xX+ x X+ yi>- y) ( B -14) 
= ^4 E <3k h «!> "IT* (2 = = . (B.15) 

k,p,q,l 

We now define a susceptibility as in the scalar case: 

XF = E Ax-y\connected (B.16) 
x 

= ^2 E ^ ~ Ak) i 2 ~ V) «»k /3k fl, a q ) - (a k a q )) , (B.17) 

k,q 

where we subtracted the disconnected component, which should vanish anyway in the symmetric 
phase (we have checked that it, indeed, does vanish). Consider also the "Majorana" correlator: 

K-y = {^ T _xl2X+xX+yl24? y ) (B.18) 

= -4^4 E (° k ft> ^ > e<<Pk W ff +P)l w iv (q+1)y (( 2 - A q)( 2 - Ai)e-^ q - AqAie"^ 1 ) , 

k,p,q,l 

and the corresponding "Majorana susceptibility:" 

XF = E ^a;-« I connected (B.19) 

= -T^E^^^-^q + ^q) (K^-k^q«-q)-(«k/3-k)(/S q a- q )) , 
k,q 

where, again, we subtracted the disconnected part, which (as we checked, once more) vanishes in 
the symmetric phase. 
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N h k Metropolis Wolff + det RW 



4 4 0.1 1.423(5) 1.424(8) 
4 2 0.5 3.65(4) 3.57(4) 
8 4 0.1 1.026(4) 1.021(8) 

Table 2: Sample comparison points, Metropolis results versus Wolff algorithm with determinant 
reweighting. 



C Simulation details 

The configurations of XY fields 77 are generated using the Wolff single-cluster algorithm [40] , 
which is well-known to overcome critical slowing down for this lattice system. The fermion measure 
is taken into account through determinant reweighting: 

(Odet M) n 

Here, O is any observable, and (• • •) denotes an expectation value with respect to the measure of 
the XY model (cf. eq. (PZTTT]) with U = 1), 

dn(rj) = Z x l Y (j\ dr,}j exp(-S K ). (C.2) 

We monitor the reliability of this method in three ways. 

1. We measure the autocorrelation time for reweighted quantities (O det M) , as well as (det M) , 
to be certain that the configurations remain independent with respect to the new measure. 

2. We perform a jackknife error analysis of the averages (O) that are obtained, gathering suffi- 
cient data to keep errors small. (An "overlap problem" would be indicated by large errors.) 

3. We have simulated a number of sample points in parameter space by an alternative method 
that does not rely on determinant reweighting. Namely, we have used Metropolis updates 
that include AS = — In det M in the action. We check that this alternative (significantly 
slower, due to much longer autocorrelation times that occur when using the Metropolis algo- 
rithm) method leads to results that agree with the determinant reweighting method. Sample 
comparison points are presented in Table [2j 

We find that the Wolff cluster algorithm together with determinant reweighting is efficient and 
reliable for all quantities and regions of parameter space that we have explored. The only difficulty 
that occurs is a sign problem in the h < 1 regime, as mentioned in the main text. However, this is 
easy to detect due to the associated large statistical errors. 

C.l Comments on y~ l 7^ 

To include y < 00 in our simulations, we have to work with a Pfaffian, as is clear from the 
results of Section (13. 2|) . eqns. (13.171 l332j) . A method to compute the Pfaffian, including the phase, 



29 



is found in the appendix of ref. [42]. We want to determine the sign ambiguity that occurs in 
Pf A4 = idet.A/1 1 / 2 , since it is obviously crucial to any averaging; here, we denoted by A4 the 
fermion matrix determined by eqns. (|3.17l [3722]) at y < oo. Similarly, we denote by M the fermion 
matrix in the y — > oo limit. 

For jT 1 < 1CT 3 , we find that InPf M = lndet M to within 4 or 5 digits for N = 4 and N = 8 
lattices. That is so close that it will not change any of the results that were stated in Sections (13.31 
13. 4p above. The result In Pf M. = In det M does not appear to depend on the lattice size. 

Right at h = 1, M becomes singular for arbitrary <p backgrounds and there is an increased 
sensitivity to the y" 1 corrections. 

Returning to h ^ 1, here are some more empirical results. For y _1 <; 0.1, Pf M is noticably 
different from det M. However, the complex phase seems much less sensitive to y^ 1 corrections than 
the magnitude. The complex phase of Pf A4 does fluctuate a bit from configuration to configuration. 
But for y = 10, it is a fluctuation in the fourth significant digit — hardly a "complex phase problem." 
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Figure 1: Comparison of susceptibilities for n = 0.1. The dashed line indicates the susceptibilities for the 
pure XY model with same k (undistinguishable, within errors for N = 4,8, 16). Here and below, large errors 
at h = 0.7 and 0.8 are due to the sign problem at h < 1. 
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Figure 2: Comparison of susceptibilities for n = 0.5. The dashed lines indicate the susceptibilities for the 
pure XY model with same k. 
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Figure 3: Comparison of susceptibilities for n = 0. The results are qualitatively similar to those obtained 
for k = 0.1,0.5. 
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Figure 4: Comparison of Binder cumulant for n = 0.1. The dashed lines indicate the values of the Binder 
cumulant for the pure XY model with same k. 
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Figure 5: Comparison of Binder cumulant for n = 0.5. The dashed lines indicate the values of the Binder 
cumulant for the pure XY model with the same n. 




Figure 6: The vortex density for k = 0.5 as a function of h. The pure XY model values of the vortex 
densities for k = 0.5 are shown in the corresponding dashed lines. The slight increase of the vortex density 
at h > 1 (especially for N = 8, 16) implies that the fermions lower the vortex energy. 
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Figure 7: The magnitude of the "Dirac" susceptibility xf for k = 0.5. The values of all physical sus- 
ceptibilities are obtained from the plot by multiplying the plotted value by the inverse square power of the 
dimcnsionless Yukawa coupling 1/Y 2 = l/(ay) 2 . 
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Figure 8: The magnitude of the "Dirac" susceptibility xf for k = 0.1. 
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Figure 9: The argument of the "Dirac" susceptibility xf for n — 0.5. 
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Figure 10: The argument of the "Dirac" susceptibility xf for k = 0.1. 
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Figure 11: The magnitude of the "Majorana" susceptibility \' F for k = 0.1. 
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Figure 12: The magnitude of the "Majorana" susceptibility x'f f° r K = 0-5. 
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Figure 13: The argument of the "Majorana" susceptibility x'f f° r K = 0.5. 
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Figure 14: The argument of the "Majorana" susceptibility x'f for re = 0.1. 
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Figure 15: The lower bound on the charged mirror inverse propagator eigenvalue, ./Vmin Q 1 , for k = 0.5. 
The dashed lines show the same quantities in the broken (n — ► oo) phase; their oscillations and dips are 
explained by eqn. (|B.10|) . The minimum mass eigenvalues are obtained, in units of L~ hys , by multiplying 
the plotted quantity by the dimensionless (large) Yukawa coupling Y\ see Section l3~4l 
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Figure 16: The lower bound on the neutral mirror fermion eigenvalue, ./Vmin f2 2 , for k = 0.5. Note that 
the charged fermions are heavier than the neutral in the symmetric phase. 
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Figure 17: The lower bound on the charged mirror fermion eigenvalue, TVmin CI 1 , for k = 0.1. 
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Figure 18: The lower bound on the neutral mirror fermion eigenvalue, TVmin fl 2 for k = 0.1. 
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